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Abstract 

We study the noise effects in a driven system of globally coupled oscillators, 
with particular attention to the interplay between driving and noise. The 
self-consistency equation for the order parameter, which measures the collec- 
tive synchronization of the system, is derived; it is found that the total order 
parameter decreases monotonically with noise, indicating overall suppression 
of synchronization. Still, for large coupling strengths, there exists an opti- 
mal noise level at which the periodic (ac) component of the order parameter 
reaches its maximum. The response of the phase velocity is also examined 
and found to display resonance behavior. 
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I. INTRODUCTION 



The set of coupled nonlinear oscillators serves as a prototype model for a variety of 
self-organizing systems in physics and in other sciences, which display the remarkable phe- 
nomena of collective synchronization [0-^. Due to analytic simplicity and some physical 
as well as biological apphcations, the system with global coupling has been mostly stud- 
ied, both analytically and numerically f^-^. Here external periodic driving may induce 
characteristic mode locking of each oscillator, leading the system to display periodic syn- 
chronization |TD|. In such a driven system, the presence of noise raises another interesting 
possibility of stochastic resonance (SR), which lead to amplification of the response of the 



system by cooperative interactions between the noise and external periodic driving ||lT|. The 
SR phenomena, which have various practical applications [p!^ -[T7[|, have been investigated in 
systems with relatively few degrees of freedom, and observed in bistable systems and also 
in systems with periodic potentials |]I4[. On the other hand, the SR effects have hardly 
been examined in a system with many degrees of freedom such as the system of coupled 



oscillators 18 



In this paper we consider a system of globally coupled stochastic oscillators, driven peri- 
odically, and investigate the interplay of noise and periodic driving, with particular attention 
to the possibility of stochastic resonance. For this purpose, it is crucial to consider appro- 
priate responses of the system to the periodic forcing. Here we consider the response of the 
phase velocity as well as the order parameter which describes the phase synchronization. We 
first derive the self-consistency equation for the order parameter and investigate the behavior 
of the order parameter in the presence of noise. It is found that the total order parame- 
ter, which consists of the time-independent (dc) and periodic (ac) components, decreases 
monotonically with noise, indicating the overall suppression of phase synchronization. The 
ac component, on the other hand, may first increase as noise grows from zero, and reach its 
maximum at a finite noise level. Such SR-like behavior is also observed in the response of 
the phase velocity: At low noise levels, the noise subtracted power spectrum of the phase 
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velocity tends to increase with noise. 

This paper consists of six sections: Section II introduces the driven system of coupled 
oscillators subject to random noise. The recurrence relation for the Fourier components is 
obtained. In Sec. Ill, we use the recurrence relation obtained in Sec. II, and derive the self- 
consistency equation for the order parameter. The set of coupled equations of motion for 
the system is transformed into a Fokker-Planck equation and the corresponding probability 
density is expanded as a Fourier series. Sections IV and V are devoted to the investigation of 
the responses of the phase and of the phase velocity, respectively, to the external driving. In 
spite of the overall suppression of synchronization, the ac component of the order parameter, 
corresponding to the phase response, as well as the response of the phase velocity is revealed 
to display SR-like behavior. Finally, a brief summary is given in Sec. VI. 

II. DRIVEN SYSTEM OF COUPLED OSCILLATORS 

The set of equations of motion governing the dynamics of the system of N coupled 
oscillators is given by 
. K ^ 

<^i + ]yEsin(0i-0i) =^i + ^iCOsQi + ri(t) {i = l,2,...,N), (1) 

where 0, represents the phase of the ith oscillator. The second term on the left-hand side 
corresponds to the global coupling between oscillators, with strength K/N . The first and 
the second terms on the right-hand side describe the natural frequency of the ith oscillator 
and the periodic driving on the zth oscillator, respectively. Finally, Ti{t) is independent 
white noise with zero mean and correlation 

{Ti{t)T^{t')) = 2D5,i5{t-t'), (2) 

where D {> 0) plays the role of the "effective temperature" of the system. The natural 
frequency uji is distributed over the whole oscillators according to the distribution g{uj)i 
which is assumed to be smooth and symmetric about ujq. Without loss of generality, we may 
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take ujq to be zero and assume that g{uj) is concave at u = 0, i.e., g"{0) < 0. The periodic 
(ac) driving amphtude /j may also vary for different oscillators, while the frequency Q of 
the driving is assumed to be uniform for all oscillators. In the absence of noise {D = 0), 
Eq. (|I|) precisely reduces to the set of equations of motion studied in Ref. [|T0|. The set 
of equations of motion in Eq. (|I]) describes a superconducting wire network and may 
also be regarded as the mean-field version of an array of resistively shunted junctions, which 
serves as a common model for describing the dynamics of superconducting arrays [Q. In 
these cases, the two terms on the right-hand side of Eq. (|l|) correspond to the combined 
direct and alternating current bias. 

Collective behavior of such an A^-oscillator system is conveniently described by the com- 
plex order parameter 

1 ^ 

^ = ]vEe^^^ =^e''' (3) 

where nonvanishing \I' indicates emergence of synchronization. Note that the synchronized 
state corresponds to the superconducting state with global phase coherence in the case of 
the superconducting network or array. The order parameter defined in Eq. (^) allows us to 
reduce Eq. into a single decoupled equation 

(pi + KA sin(0j — 6) = uJi + li cos Qt + Ti{t). 

We then seek the stationary solution with 6 being constant, which is possible due to the 
symmetry of the distribution of Ui and li about zero. Redefining (pi— 6 as (pi and suppressing 
indices, we obtain the reduced equation of motion 

(p + KAsm(p = uj + Icosnt + T{t), (4) 

which depends explicitly on the order parameter. In this manner the order parameter A, 
defined in terms of the phase via Eq. (^, in turn determines the behavior of the phase via 
Eq. (^), and can thus be obtained by imposing self-consistency, as discussed in Sec. III. Note 
here that A in general depends periodically on time due to the periodic driving; this allows 
the Fourier expansion 



A = Ao + A, cos{snt + as) (5) 

s=l 

with appropriate phases as, where Aq is the time- independent (dc) component and A^ is 
the time-dependent (ac) one due to the periodic forcing. 

A convenient way to deal with a set of Langevin equations is to introduce an appropriate 



probabihty density and to resort to the associated Fokker-Planck equation In general 
the set of Langevin equations (|l]) makes it necessary to consider the A^-oscillator prob- 
ability density P{{(j)i},t) and the corresponding Fokker-Planck equation 0. In the system 
with glabal coupling, however, the set in Eq. (0) naturally reduces to the single Langevin 
equation (^, as shown above. This in turn leads to the Fokker-Planck equation for the 
single-oscillator probability density P{4>, t) with the self- consistency for the order parameter 
explicitly imposed, which has been considered in the absence of driving [§]. 



The Fokker-Planck equation for the probability density P{4>, t) reads ^1 



dP _ d 
'dt ~dd> 



d^P 

D^TT, (6) 



where V{(j)) = —KAcos(j) — is the washboard potential. Unlike the system without 
driving (/ = 0), the stationary solution of which has been obtained [Q, Eq. does not 
allow such a simple stationary solution. We thus use the periodicity of the system and 
expand the probability density Fourier series 

oo 

P(0,t)= C^n(^)e^"^ (7) 

n=— oo 

which, upon substitution into Eq. (||), yields 

Cn{t)= - [in{uj + / cos fit) + u'^d] Cnif) 

- -KA{t)Cn+i{t) + -irA(t)a-i(t). (8) 

Since the probability density should be real, we have the relation C„ = Cl„; the normal- 
ization condition Jq^ P{(f),t)d(j) = 1 gives the constant term Co = l/27r. The differential 
recurrence relation in Eq. (Rl) can be written in the form of an integral recurrence equation 



C„(t)= CniO) exp 



-K exp 



-in{ijjt + — sin Vtt) — n^Dt 
—in{ujt + — sin Vtt) — n^Dt 



X / dt'/\{t') Cn+l{t') - Cn-l{t') 

Jo 



exp 



in{ujt' 



n 



sinQt') + n^Dt' 



(9) 



which is of the same form as the equation for a single oscillator p2|, except for that here 
self-consistency for the order parameter is required. 



III. SELF-CONSISTENCY EQUATION FOR THE ORDER PARAMETER 

In this section we derive the self-consistency equation for the order parameter, which 
describes the response of the phase and determines the collective behavior of the system. We 
suppose that the periodic driving amplitude / is distributed according to /(/), independently 
of the natural frequency u. Recalling that in Eq. in fact represents (f) — 0, we have the 
self-consistency equation 

N ^ 

3 

/OO /"OO 
/ dujg{uj){e^^)^^i, (10) 
-OO J —OO 

where (■ ■ denotes the average with given uj and I. 

With the probability density P{(f),t), the expansion of which is given by Eq. (|^, we 
compute the average 

/•2-K 

{e'^)^ d<Pe'^P{<P,t) = 27rCl{t), (11) 
Jo 

where the relation C„ = C*„ has been used. Then Eq. (p!0| ) leads straightforwardly to 

dlfil) / du;giu)Clit), (12) 

-OO J —OO 

which gives the order parameter in terms of the Fourier coefficient Ci [t) . Assuming KA ^ 1 
near the transition to the coherent state, we need to obtain Ci up to the order of (KA)^. 
For this, we first compute C2 from Eq. (||), neglecting C3, and substitute the obtained C2 
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back into the equation for Ci, i.e., Eq. (|]) with n = 1. At long times the transient terms 
such as exp[—n'^Dt] for n 7^ vanish, and a lengthy calculation yields 



■iK °° 



Qi{£—m+s)Qt+ias ^i{l—Tn—s)Qt—ias 



iu+{£ + s)n-iD iu + {£-s)n-iD 

~ 51 As^s'^s" Yl Jl{x)J(,i{x)Jm{2x)Jm'{2x)Jn{x)Jn'{x) 

s,s',s" £,£' ,m,,Tn' ,n,n' 

X [e*("=-"^'-"=")F(t; s, s', s") + e^("-+"-'-"-")F(t; s, -s', s") 
+ e^("--"s'+"«")F(t; -s, s', s) + e*(°=+"='+°^")F(t; -s, -s', s) 
+ e^(-"^-"«'-"=")F(t; s, s', -s) + e^(-°^+".'-°.")F(t; s, -s', -s) 
+e*(-"^-"^'+"^")F(t; -s, s', -s") + e^(-"-°s'-°s")F(t; -s, -s', -s")] , (13) 

where x = I /VL, = 0, and the function F{t] s, s', s") depends on the indices £, i', m, m', n, 
and n' as well as on u and Q: 

+m-m' +n-n' -s-s' ~s")Q.t 

F{t;s,s', s": 



[iu + {n' + s)il - iD] [2iu + {m' - n + n' + s + s')il - AiD] 

1 



u! + {£' — m + m' — n + n' + s + s' + s")^l — iD 

With the above expression for Ci, Eq. (|12D yields the explicit form of the self-consistency 
equation for the order parameter. 

Comparing term by term in the resulting self-consistency equation, we can determine 
each component of the order parameter. Namely, the dc component Aq is given by the 
constant (zero- frequency) terms in the expansion of Ci{t). The next component Ai can be 
obtained from the terms with frequency Q, the component A2 from the 2Q terms, and so 
on. For weak driving, the ac components of the order parameter are much smaller than the 
dc component, leading to the simple self-consistency equation: 

A ^ aKAo - b{KAof (14) 

with the coefficients 



«= ^E/ dl f{I)Ux)Jm{x)e'^'-'^^''' j duj 

£,Tn 



uj + mVt + iD 



I 

1,1' ,m,m' ,n,n' 



b=-- E / dIf{I)Mx)Jeix)JU2x)J^i2x)Ux)Jnix) -^+n' -n)nt 



I 1 9{^> 



[uj + n'Vt + iD] [2uj + {m' + n' - n)Vt + AiD] 



X \ (15) 

u + {I' + m' - m + n' - n)VL + iD' ^ ' 

In the simple case of no external driving (/ = 0) and noise {D 0), the representa- 
tion Ti5{uj) = Diuj"^ + D^)^^ in the limit D ^ together with the symmetry of g{uj) 
reduces Eq. (|15|) to a = (7r/2)(7(0) and b = — (7r/16)(7"(0), which indeed reproduces the 
self-consistency equation obtained in Ref. 0]. 



Solving Eq. (|T^), we obtain the collective behavior of the system, which has been analyzed 
in Ref. [jlO[: For small K, only the trivial solution A = exists. On the other hand, for 
K > Kc = l/ttQ, Eq. ( ]T^ ) also allows the nontrivial solution with the dc component 



JboK{aoK - 1) 



where the constant coefficients are given by 



duo 



i,m,m'n,n' 



/ dl f{I)Ji{x)Jm{2x)Jrr,f{2x)Jn{x)Jn'{x)Ji+rn+n-m'~n'{x) 



^ I [uj + n'n + iD] [u + m + iD] [2uj + (m' + n' - n)n + AiD] ' ^^^^ 

Thus as K is increased beyond K^, the null solution becomes unstable and the stable non- 
trivial solution A+ (accompanied by the ac components) appears via a pitchfork bifurcation 



at K = Kr 10 . 



IV. NOISE EFFECTS ON SYNCHRONIZATION 

To understand the cooperative effects of the driving and noise on the response of the 
system, we examine in this section how the noise affects synchronization behavior of the 
system. For simplicity, we consider the weak-driving or high-frequency limit {x = I /Q <^ 1), 
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expand the Bessel functions in Eq. (^) to the order of x^, and perform the average over 
the distribution /(/). This gives the coefficient a to the order of ct/, the variance of the 
distribution /(/): 

D(cos 2Vtt -l) + {uo + Vt) sin 2fit L>(cos 2^]t -l)-{uj -Vt) sin 2nt 

{uj + + {uj - ny + 

^ D cos 2(]t + (cu + 21]) sin 2nt ^ D cos 2fit - - 2fi) sin 2fit 



2(cj + 2fi)2 + 2^2 2(cu - 2fi)2 + 2^2 

= ao + 02 cos(2r2t + a2)- (18) 

Similarly, a tedious but straightforward calculation leads to the coefficient b = Bq + 
62 cos(2i7t + 02). Note that the symmetry of the distribution /(J) about / = forbids 
the frequency Q term, which is linear in the driving. Here it is easy to observe that ao, 
which reads 



2 V 2^2; y Cj2 + D2 

1 



+ 



(19) 



in general decreases monotonically with D. Thus the critical coupling strength grows 
as the noise level is raised. Figure 1 displays the monotonic increase of Kc (= ao^^) with 
the noise level D, for o"/ = 0.1 and Vt = 2.0. For the distribution of natural frequencies, the 
Gaussian distribution with variance a^^ = 0.5 has been chosen. 

With the coefficients a and b obtained above, the order parameter can be obtained 
from Eq. (|l^ and its behavior in the presence of noise can be investigated. Indeed the 
dc component Ao given by Eq. (|16D is easily found to decrease monotonically as the noise 
level D is increased. On the other hand, it is too complicated to obtain analytically the 
explicit behavior of the ac component (s > 1). Further, the analytical results are based 
on Eq. (p^, which is valid only near the transition [K ^ Kc); this makes it desirable to 
obtain the order parameter numerically. We have thus performed numerical simulations to 
compute the components Aq and A2. Since the effects of external driving first appear in 
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the coefficients 02 and 62, giving rise to the frequency 2Q term, it is relevant to investigate 
A2 as the response to the external driving. In the simulations Eq. has been integrated 
with discrete time steps of At = 0.01. At each run, we have used Nt = 6048 time steps to 
compute the order parameter, discarding the data from the first 4 x 10^ steps, and varied 
both At and A^^ to verify that the stationary state was achieved. Finally, independent 
runs with 30 different distributions of the natural frequency and initial conditions have been 
performed, over which the averages have been taken. For both the distribution of the driving 
amplitudes and that of the natural frequencies, we have chosen Gaussian distributions with 
various values of variances aj and o"^, only to find no qualitative difference. 

Figure 2 displays the obtained behaviors (a) of the dc component Aq and (b) of the ac 
component A2 in the system of = 1000 oscillators, driven by frequency Q = 1.0738 and 
with variances = 1.0 and aj = 1.0. The data represented by empty and solid squares in 
Fig. 2 correspond to the coupling strength K = 2.5 and K = 3.0, respectively. It is shown 
that for both values of the coupling strength, Aq decreases monotonically as the noise level 
D is raised. Such monotonic behavior is also exhibited by A2 for K = 2.5. For K = 3.0, on 
the other hand. Fig. 2(b) displays that the ac component A2 first increases with noise and 
reach its maximum at a finite value of the noise level D. Similar non- monotonic behavior 
of A2 can be observed for larger values of the coupling strength K, suggesting the presence 
of SR-like behavior in the order parameter. Note, however, that the dc component Aq is 
in general dominant over the ac component, leading to the monotonic decrease of the total 
order parameter. It is thus concluded that noise tends to suppress monotonically the overall 
synchronization in the system. 

V. RESPONSE OF THE PHASE VELOCITY 

In this section we investigate the power spectrum of the phase velocity at the driving 
frequency, which conveniently describes the response of the phase velocity to the external 
driving. In the case of a superconducting wire network or array, the phase velocity can be 
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identified with the voltage via the Josephson relation, and the power spectrum of the phase 
velocity simply corresponds to the voltage power spectrum under the combined direct and 
alternating current driving. Equations and give the average phase velocity of a single 
oscillator in terms of ImCi, the imaginary part of Ci. 



2lT 



Jo 

= a; + / COS + 2ttKA ImCi, 



which, upon substitution of Eq. (|1^) for ImCi, obtains the simple form 



(20) 



{^) =uj + A cos nt + B sin nt + 0{K'^ AqA^ 



(21) 



with the amplitudes 



A 



B 



2VL 



uo + VL 



oj — VL 



2n 



DI 



1 2 1 



+ 



The desired power spectrum S of the phase velocity at the driving frequency is propor- 
tional to the square of the Fourier component of frequency Q, i.e., S{Q) oc A"^ + B"^. In the 
limit D ^ 0, the amplitude of the Fourier component approaches 

^TT^K^Af 



1 + 



fi2 



+ 1' 



^5{uj + fi) - 25{u}) + 5{uj - n)f 



while it approaches P in the limit D ^ oo. It is of interest to note that the amplitude in the 
noiseless limit can be either larger or smaller than that in the strong- noise limit, depending 
on oj and Vt: As the noise level is raised, the amplitude tends to decrease from the noiseless 
value for uj > Vt and increase for < Vt. Accordingly, for given driving frequency, those 
oscillators with smaller /larger natural frequencies contribute to the increase/decrease of 
the amplitude toward its asymptotic value P. For small values of the variance o"^;, for 
example, most oscillators should possess the natural frequency oj < VL, although there may 
still exist some oscillators with frequency uj > Vt. The power spectrum of the whole system, 
which is given by the sum of contributions from all the oscillators, is then expected to 
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increase for small D and to approach the asymptotic value which is proportional to P. 
Unfortunately, however, the approximations used on various stages of the analysis disallows 
a reliable analysis. In particular the extrapolation to the limit D is untrustworthy since 
nonzero effective temperature {D ^ 0) has been assumed in solving the equation for P(0, t). 
It is also obvious that the higher-order terms neglected in the analysis set a limit in the 
regime of validity, making it desirable to investigate the system by other means. 

We have thus performed numerical simulations to obtain the power spectrum for vari- 
ous values of the coupling strength and of the variance in the distributions of the natural 
frequency and of the driving amplitude. We have again integrated Eq. (|l]) for the system 
of = 1000 oscillators with discrete time steps of At = 0.01, using at each run Nt = 6048 
time steps to compute the power spectrum of the phase velocity and discarding the data 
from the first 4 x 10^ steps. The averages have been taken over 300 independent runs with 
different distributions of the natural frequency and initial conditions. /^From the obtained 
time series, we have computed the power spectrum by means of the fast Fourier transform 
algorithm. To take into account the background noise, we have taken five nearest data 
points around the peak at the driving frequency in the power spectrum and performed the 
average to give the noise level. (The results have been found not to change qualitatively 
even if other measure for the noise level is adopted.) 

In Fig. 3 we present the obtained data: the background noise subtracted power spectrum 
of the phase velocity at the driving frequency versus the noise level. For the distributions 
of the natural frequency and of the driving amplitude, Gaussian ones have been chosen 
with variances a^^ = 0.5 and aj = 0.2, respectively, while the coupling strength K = 2.5 
and the driving frequency Q = 7r/1.024 have been taken. Remarkably, Fig. 3 displays that 
the power spectrum increases as the noise level is raised from zero. Obviously it does not 
keep increasing monotonically with the noise, and there apparently exists an optimal noise 
level at which the power spectrum reaches its maximum. Beyond the optimal noise level, 
the power spectrum first falls off gradually and saturates eventually toward its asymptotic 
value, although this behavior is somewhat obscured by the large fluctuations due to strong 
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noise. Such a broad peak followed by gradual decrease has also been observed in the SR 
of another system [|I5|. It is thus suggested that the response of the phase velocity also 
displays SR-like behavior in the appropriate regime. Since the phase velocity corresponds 
to the voltage in a superconducting system, this indicates that the noise subtracted power 
spectrum of the voltage displays such resonance behavior. 

VI. SUMMARY 

We have studied the noise effects in a driven system of globally coupled oscillators, with 
emphasis on the interplay of noise and periodic driving. In particular, to investigate the 
possibility of resonance behavior, we have considered the response of the phase velocity as 
well as the order parameter which describes the phase synchronization. The self-consistency 
equation for the order parameter, derived from the recurrence relation of the probability 
density, has been shown to display monotonic decrease of the total order parameter in the 
presence of noise. It has thus been concluded that noise in general suppresses overall phase 
synchronization in the system, i.e., superconductivity tends to be disturbed by noise present 
in a superconducting wire network or array. 

Nevertheless it has also been revealed that for large coupling strengths the ac component 
of the order parameter increases with the noise level growing from zero and reach its maxi- 
mum at a finite noise level. Such resonance behavior has also been observed in the response 
of the phase velocity: At low noise levels, the noise subtracted power spectrum of the phase 
velocity has been found to increase with noise, displaying a broad peak at a finite noise 
level. As the noise level is raised further, the power spectrum appears to saturate toward its 
asymptotic value, although concealed by large fiuctuations due to strong noise. In conclu- 
sion, the phase synchronization, describing the collective behavior of the coupled-oscillator 
system, is suppressed monotonically in the presence of noise. Still the responses of the phase 
and of the phase velocity can display non-monotonic resonance behavior in the appropriate 
regime, which may be manifested by a broad resonance peak of the voltage power spectrum 



13 



in the case of a superconducting system. 

It is also of interest to note that the phase velocity on average may serve as a measure 
of A/ A, the rate of change of phase synchronization. Accordingly, the resonance behavior 
in the response of the phase velocity suggests that the approach to the coherent state with 
synchronization can be accelerated by the presence of weak noise. Since the phase synchro- 
nization corresponds to the memory retrieval in the network of neuronal oscillators 
the resonance behavior may also imply the information processing assisted by weak noise 
in a biological system, e.g., the crayfish who appears to use such resonance to perceive an 
enemy quickly. The detailed understanding of the resonance behavior and its implications to 
applicable physical and biological systems require more extensive analytical and numerical 
investigations, which are left for further study. 
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FIGURES 




FIG. 1. The critical coupling strength beyond which synchronization appears versus the noise 
level in the system with the driving frequency Q, = 2 and the variances a^, = 0-5 and aj = 0.1. The 
random noise in the system monotonically increases the critical coupling strength, thus tending to 
suppress synchronization. 
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FIG. 2. Behavior of the order parameter in the presence of noise. The data represented by 
empty and sohd squares correspond to the couphng strength K = 2.5 and K = 3.0, respectively, 
(a) The dc component Aq for both K = 2.5 and K = 3.0 is shown to decrease monotonically. (b) 
The ac component A2 for K = 2.5 decreases monotonically, while for K = 3.0 it displays a peak at 
a finite noise level. The standard deviations of the data (not shown) range from 5% to 15%, and 
lines are merely guides to the eye. 
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FIG. 3. The noise subtracted power spectrum of the phase velocity at the driving frequency. 
There appears an optimal noise level at which the power spectrum reaches its maximum. The error 
bars have been estimated by the standard deviation and the line is merely a guide to the eye. 
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